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ABSTRACT 

The  development  of  fully  nonlinear  Gortler  vortices  in  the  high  Reynolds  number  flow 
in  a  symmetrically  constricted  channel  is  investigated.  Attention  is  restricted  to  the  case  of 
‘strongly’  constricted  channels  considered  by  Smith  and  Daniels  (1981)  for  which  the  scaled 
constriction  height  is  asymptotically  large.  Such  flows  are  known  to  develop  a  Goldstein 
singularity  and  subsequently  become  separated  at  some  downstream  station  past  the  point 
of  maximum  channel  constriction.  It  is  shown  that  these  flows  can  support  fully  nonlinear 
Gortler  vortices,  of  the  form  elucidated  by  Hall  and  Lakin  (1988),  for  constrictions  which 
have  an  appreciable  region  of  local  concave  curvature  upstream  of  the  position  at  which 
separation  occurs.  The  effect  on  the  onset  of  separation  due  to  the  nonlinear  Gortler  modes 
will  be  discussed.  A  brief  discussion  of  other  possible  nonlinear  states  which  may  also  have 
a  dramatic  effect  in  delaying  (or  promoting)  separation  will  be  given. 
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DTIC  QUALITY  INSPECTED  1 


1.  Introduction 


Our  concern  is  with  the  development  of  fully  nonlinear  Gortler  vortices  in  Smith- 
Daniels  flows  and  their  effect  on  the  onset  of  separation.  By  Smith-Daniels  flows  we 
refer  to  the  steady  flow  in  a  symmetrically  constricted  channel  considered  by  Smith  and 
Daniels  (1981)  who  showed  that  when  the  scaled  height,  h  1,  of  the  constriction  becomes 
large  a  secondary  classical  boundary  layer  of  thickness  0{h~1!2)  is  set  up  within  the  wall 
boundary  layer  (of  thickness  0(Re-1/3)  where  Re  is  the  global  Reynolds  number  for  the 
flow).  This  inner,  classical  boundary  layer  develops  a  Goldstein  (1948)  singularity  at  a 
position  downstream  of  the  position  of  maximum  channel  constriction.  However,  Smith 
and  Daniels  demonstrated  that  this  singularity  could  be  removed,  without  any  upstream 
influence  being  induced,  and  their  analysis  was  able  to  describe  the  reattachment  of  the 
boundary  layer  at  a  distance  0(h 3)  downstream. 

Subsequently  Hall  and  Bennett  (1986)  considered  the  linear  instability  of  this  flow, 
in  the  general  context  of  the  stability  of  both  steady  and  unsteady  interactive  boundary 
layer  flows,  at  stream-wise  locations  before  the  Goldstein  singularity  develops.  This  study 
demonstrated  that  the  flow  is  unstable  to  Gortler  vortices  for  channels  in  which  the  con¬ 
striction  has  a  region  of  concave  curvature,  in  an  interval  before  the  position  at  which 
separation  occurs.  In  this  analysis  fxx,  where  f(x)  =  hF{x)  with  h  >  1  is  the  constriction 
profile,  effectively  plays  the  role  of  the  Gortler  number. 

In  the  case  of  external  boundary  layer  flows  over  a  wall  of  variable  curvature  Hall 
and  Lakin  (1988),  motivated  by  the  earlier  weakly  nonlinear  theory  of  Gortler  vortices  by 
Hall  (1982b)  and  the  numerical  solution  of  the  fully  nonlinear  Gortler  vortex  governing 
equations  by  Hall  (1988),  developed  a  self  consistent  asymptotic  theory  for  short  wave¬ 
length,  large  amplitude  Gortler  vortices.  In  this  theory  the  Gortler  number  G  scales  with 
k 4  (where  A:  1  is  the  vortex  wave  number).  Their  analysis  demonstrated  that  the  mean 
flow  adjusts  due  to  the  presence  of  the  vortex  so  as  to  render  all  vortex  modes  neutrally 
stable.  Within  an  0(1)  region  of  the  boundary  layer  the  mean  flow  is  driven  by  the  vortex 
velocity  field  and  as  such  has  no  relation  to  the  unperturbed  boundary  layer  flow.  Such  a 
situation  is  similar  to  that  postulated  by  Malkus  (1956)  for  turbulent  flows;  Malkus  argued 
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that  the  mean  part  of  a  turbulent  flow  would  adjust  so  as  to  ensure  that  all  ‘modes’  are 
rendered  neutrally  stable. 

As  pointed  out  by  Hall  and  Lakin  (1988),  the  modification  of  the  mean  flow  by  the 
nonlinear  longitudinal  vortex  motion  raises  the  question  as  to  what  effect  this  modified 
mean  flow  has  on  the  onset  of  separation,  in  separating  flows  that  can  support  Gortler 
vortex  like  disturbances.-  One  important  result  of  the  Hall  and  Lakin  calculation  is  that 
the  skin  friction  for  the  vortex  induced  mean  flow  can  be  dramatically  increased  over  that 
of  the  unperturbed  boundary  layer  (although  no  quantitative  results  on  this  point  are 
presented  this  can  be  inferred  from  their  figure  6,  Hall  and  Lakin  (1988)  pp  440-441);  we 
would  then  expect  that  the  onset  of  separation  could  be  considerably  delayed  due  to  the 
presence  of  the  vortex  motion.  In  light  of  the  work  of  Hall  and  Bennett  (1986)  alluded  to 
above,  we  choose  to  concentrate  on  the  Smith  and  Daniels  (1981)  problem  of  the  steady 
flow  in  a  symmetrically  constricted  channel. 

The  procedure  adopted  for  the  rest  of  this  paper  is  as  follows.  In  §2  we  develop 
the  governing  equations  for  nonlinear  Gortler  vortices  in  the  0(h~ 1/f2)  classical  boundary 
layer  of  Smith  and  Daniels  together  with  a  discussion  of  the  relevant  results  from  Smith 
(1982)  concerning  the  onset  of  separation  as  well  as  the  those  points  from  the  linear  theory 
of  Hall  and  Bennett  (1986)  pertinent  to  the  present  investigation.  In  §3  we  extend  the 
linear  theory  of  Hall  and  Bennett  (1986)  to  the  fully  nonlinear  regime  for  short  wavelength 
vortices  considered  by  Hall  and  Lakin  (1988).  In  §4  we  present  the  numerical  solution 
of  the  vortex  induced  mean  flow  equations.  Finally  in  §5  a  discussion  of  our  results  is 
presented;  particular  emphasis  is  given  to  the  effect  that  the  presence  of  the  nonlinear 
vortex  structure  has  upon  the  onset  of  separation. 

2.  Formulation  of  the  governing  equations 

Consider  the  steady  motion  of  an  incompressible  fluid  of  density  p  and  kinematic 
viscosity  u  in  a  symmetrically  constricted  channel.  Let  L  denote  the  typical  length  scale  of 
the  constriction  and  let  Uoo  denote  the  typical  free-stream  speed  sufficiently  far  upstream 
of  the  constriction  (see  figure  1).  We  define  the  Reynolds  number  Re  by  Re  =  UocL/v, 
and  throughout  this  work  we  will  assume  that  the  flow  is  such  that  the  Reynolds  number 
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is  asymptotically  large.  The  steady  three-dimensional  Navier-Stokes  equations  in  nondi- 
mensional  form  are 


du  dv  dw 

dx  dy  **"  dz  ' 

du  du  du  dp  / d2u  d2u  d2u\ 

Ud^  +  Vd^  +  Wdl =  ~lh+Re  [d^2  +  d^2  +  dl2 )  ' 

dv  dv  dv  dp  / d2v  d2v  d2v\  (^-l) 

Ulh+Vfy+Wlh=:~fy+  e  \JW  +  Ihf  +  dl*  )  ' 

dw  dw  dw  dp  _  / d2w  d2w  d2w\ 

Ulh+Vlh;  +  W7h=~jh+Re  V Iti2  +  W  +  '&?)  ’ 

where  (r,  y,  2)  are  the  usual  Cartesian  co-ordinates,  nondimensionalized  with  respect  to 
X,  ( u,v,w )  are  the  corresponding  velocity  components,  nondimensionalized  with  respect 
to  the  upstream  velocity  Uoo  and  p  is  the  pressure,  nondimensionalized  with  respect  to 
pU^.  For  definiteness  we  will  assume,  following  Smith  (1982),  sufficiently  far  upstream  of 
the  channel  constriction  the  flow  is  Poiseuille 


u  -*  y  -  y2,  v  -*  0,  p  -*■  - 


x  —00. 


As  shown  by  Smith  (1982),  with  a  constriction  length  L  =  0(1),  the  critical  height  of  the 
constriction  is  0(Re~1/3)  (Smith’s  ‘strong  constriction’).  We  then  take  the  constriction 
y  =  hRe~' S3 F(x)  at  the  lower  wall  (a  corresponding  constriction  y  =  1  —  hRe  l/3F(x) 
is  also  present  at  the  upper  wall,  however  here  we  are  predominantly  concerned  with  the 
boundary  layer  flow  in  the  vicinity  of  the  lower  wall  and  as  such  we  will  restrict  our 
attention  to  the  lower  boundary);  we  will  subsequently  take  the  limit  h  >■  1  (this  being 
the  appropriate  asymptotic  regime  of  the  Smith  and  Daniels  (1981)  investigation)  however 
we  initially  consider  0  <  h  <  00.  The  flow  develops  a  viscous  wall  layer  of  thickness 
0(Re~1/3),  in  which  the  velocity  field  and  pressure  expand  as  (Smith  (1982)) 

(u,  v,  id, p)  =  (Re~'/3U,  Re~2/3V ,  Re~2'3W ,  Re~2/3P0  +  Re~4/3Px ), 

(2.2) 

x  —  0(1),  y  =  Re~l'3Y,  z  =  Re~l'3Z , 
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where  Po  is  a  function  of  x  while  P\  is  a  function  of  x,Y,Z.  The  governing  equations  in 
the  viscous  wall  layer  are  then  given  by 


dU  dV  dW  n 

dx  *  dY  *  dz 

udJL^vdJL  +  wdJL  =  JJi  +  d2lL  +  d^L 

dx  dY  dZ  dx  dY2  dZ 2  ’ 

uW  +  vTL  +  w™ 

dx*  dY  dZ  dY  dY 2  dZ 2  5 

dW  dW  dW  dPx  d2W  d2W 

U  dx  +V  dY  dZ  ~  dZ  +  dY2  +  dZ2  ’ 

together  with  the  boundary  conditions 

(u,u,  w,p)  — >  (Y,  0,0,0)  as  x  — >  — oo, 

u  =  v  =  w  =  0  on  Y  =  hF(x ), 
f/~F  F  ->  oo. 


(2.3) 


(2.4) 


Here  we  have  anticipated  the  well  known  result  that  the  Gortler  vortex  wavelength  is  scaled 
on  the  boundary  layer  thickness. 

In  the  asymptotic  limit  h  »  1,  Smith  and  Daniels  (1981)  have  demonstrated  that  the 
flow  described  above  develops  a  classical  boundary  layer  of  thickness  0(h~ l!2)  attached 
to  the  hump.  In  the  absence  of  any  vortex  component  the  appropriate  scalings  axe  found 
to  be 

U  =  hu,  V  =  h}!2v,  Po  =h2p  £  =  h}<2Y. 


From  the  analysis  of  Smith  and  Daniels  (1981)  the  governing  equations  are  given  by 

du  dv 

w;-  al  +  a«=0’  <2-5o’i,) 

together  with  the  boundary  conditions 


_du  _du  d2u 

udi  +  vdl~  I+a?’ 


u  =  u  =  0,  £  =  0,  u  —*  F  £  — >  oo, 


(2.6a,  b) 


and  the  upstream  conditions 


u,  v  — ►  0  x  — >  — oo; 


(2.6c) 


4 


(here  we  axe  assuming  F(x )  — >  Oasx  — >  — oo  and  we  have  made  use  of  the  Prandtl 
transformation).  We  note  that  in  (2.5a)  the  pressure  gradient  term  —FFX  arises  from  the 
outer  flow  solutions  (the  reader  is  referred  to  Smith  and  Daniels  (1981)  for  full  details). 
Subsequently  we  will  assume  that  the  constriction  profile  is  such  that  F(x)  =  0  for  x  <  0 
so  that  the  upstream  condition  (2.6c)  can  be  replaced  by 

u,  v  =  0  x  <  0.  (2.6 d) 


To  derive  the  governing  equations  for  the  Gortler  vortex  velocity  field,  which  we 
assume  is  confined  to  the  secondary  0(h~1^2)  boundary  layer,  the  following  scalings  are 
employed  (see  Hall  and  Bennett  (1986)).  The  x,  y  and  z  disturbance  velocity  components 
scale  with  the  x,  y  and  y  basic  state  velocity  components;  the  y  and  2  variations  are  on 
the  classical  boundary  layer  scale,  namely  0(/i-1/2),  and  the  pressure  perturbations  are 
reduced  in  size  until  the  pressure  gradient  and  viscous  terms  in  the  y  and  z  momentum 
equations  become  comparable.  Thus 


i  =  h1/2F,  2  =  hx'2Z, 

U  =  hu,  V  =h}/2v,  W  =  hl/2w,  P0  =  h2p,  Pi  =  hp. 


(2.7) 


The  governing  equations  are  then  given  by 


du  dv  dw 

dx  +  ae  +  aT  =  °’ 

du  du  du  d2u  d2u 

Ufa+Vdt+W1h  =  Fx  +  de+d^' 

dv  dv  dv  3/2  2  dp  d2v  d2v 

d2w 


dx 


dw  dw  dw  dp  d2w 
Ulh  +  V~dt  +  W~d~z  =  ~~d~z  +  ~d^+  dz 2  ’ 


together  with  the  boundary  conditions 


u  —  v  =  w  —  0  ^  =  0, 

u  — *  F,  w  — ►  0  £  — ♦  oo, 

u,v,w  —  0  x  <  0, 


(2.8a  -  d) 


(2.9  a  —  c) 


where  we  have  made  use  of  the  Prandtl  transformation. 
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Our  concern  is  with  the  solution  of  the  nonlinear  system  (2.8,9)  in  the  limit  of  short 
spanwise  wavelength,  d/dz  1.  Before  proceeding  with  a  description  of  the  solution  of 
these  equations  we  make  some  comments  concerning  the  unperturbed  system  (2.5,6)  and 
recall  the  relevant  details  from  the  linear  theory  of  Hall  and  Bennett  (1986). 

In  practice  the  solution  of  the  system  (2.5,6)  is  a  numerical  task;  we  briefly  describe 
the  numerical  technique  used  to  solve  this  system.  In  order  to  obtain  solutions  of  (2.5a,c) 
to  within  the  numerical  accuracy  required  for  the  subsequent  calculations  of  §4  it  was 
found  to  be  advantageous  to  solve  a  modified  system  obtained  from  (2.5)  by  writing 


u  =  Fu0,  v  =  -Fx(  +  v0, 


where  («o,  vo)  are  governed  by 

F.al  +/i,|  +  (So  -  =  F,  + 

F'l-~  l)+^  +  f  =°’ 

with  the  appropriate  free-stream  boundary  condition  now  given  by 


Uq 


Z  -*•  °o; 


all  other  boundary  conditions  remaining  unchanged.  This  modified  system  was  discetized 
according  to 


Funun 


1  Fu 2 

— (un+i  -  2 un  +  un- 1)  =  4-  Fx{  1  -  «i) 

(vn  F t(n)(Un+l  Un  —  j) 


2  h 


(2.10a,  b) 


^n-f  1  ~  Vn-l  _  F(un  -  Un) 
2  h  e 


-  Fx(un  -  1), 


where  a  tilde  denotes  a  quantity  evaluated  at  x  +  e  and  the  index  n  refers  to  a  quantity 
evaluated  at  the  grid  point  (  =  nh;  we  explicitly  impose  the  boundary  conditions  (2.6a,b,d). 
Given  u„,  vn  and  F  the  tridiagonal  system  (2.10a)  can  be  solved  for  the  updated  value  un; 
vn  is  then  determined  from  (2.10b).  The  solution  was  marched  downstream  from  x  =  0, 
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with  the  initial  conditions  (2.6d),  iterating  at  each  step  until  converged  values  of  u  and  v 
were  obtained.  For  definiteness  we  will  restrict  our  attention  to  the  hump  profile 

0,  x  <  0, 

F(x)  =  <  ae~1(x4/ 16  —  x3  +  Ax2),  0  <  x  <  4,  (2-11) 

ax2e-;c2/16,  x  >  4, 

where  a  is  an  amplitude  parameter.  Consider  now  the  boundary  conditions  (2.9b).  As 
noted  by  Smith  (1982)  the  boundary  layer  is  accelerated  where  Fx(x)  >  0  (upstream 
of  the  point  of  maximum  constriction  x  —  xmax  =  4).  Downstream  of  xmax  the  flow  is 
decelerated  due  to  the  adverse  pressure  gradient  —FFX\  the  skin  friction  then  tends  to  zero, 
the  Goldstein  singularity  is  encountered  and  the  boundary  layer  flow  becomes  separated. 
In  figure  2  we  present  a  plot  of  the  skin  friction  t(x)  =  ii{(x,  0)  versus  x  for  the  profile 
(2.11)  and  the  particular  choice  ft  =  1.031.  We  see  that  r  attains  a  maximum  at  x  «  2.9 
and  then  decreases  until  r  =  0  at  x  =  xs  «  4.53  at  which  point  the  Goldstein  singularity 
will  arise  and  separation  will  occur.  We  note  that  changing  the  amplitude  ft  has  no  effect 
on  either  the  position  at  which  r  attains  its  maximum  nor  the  point  at  which  r  =  0.  Our 
interest  lies  in  the  effect  that  the  presence  of  a  fully  nonlinear  vortex  structure  upstream 
of  xmax  may  have  in  delaying  (or  promoting)  the  onset  of  separation. 

The  linear  theory  of  Hall  and  Bennett  (1986)  demonstrates  that  the  boundary  layer 
becomes  unstable  to  linear  Gortler  vortices  at  the  position  x*  where 

du 

=  1,  (2.12) 

and  are  confined  to  a  thin  layer  about  the  position  £  =  £*  where 

d  / _du\ 

dtVdt) 

As  noted  by  Hall  and  Bennett,  for  humps  with  F  >  0  (ie  channel  constriction)  uu^  remains 
positive  until  the  Goldstein  singularity  develops,  at  x  =  x,  say,  and  hence  for  instability  we 
require  the  constriction  to  be  locally  concave  somewhere  in  the  interval  0  <  x  <  xs.  The 
profile  (2.11)  was  chosen  with  this  in  mind;  it  is  locally  concave  in  the  region  0  <  x  <  1.69. 
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As  an  aside,  we  note  that,  from  the  results  of  Hall  and  Bennett,  there  is  a  critical 
amplitude  a  below  which  instability  cannot  occur  (see  expression  (3.15)  and  the  preceding 
discussion  in  Hall  and  Bennett  (1986)).  In  figure  3  we  present  a  plot  of  the  neutral 
position  x*  versus  the  amplitude  a  for  the  profile  (2.11);  we  find  the  critical  amplitude  to 
be  a  =  ac  ss  0.84.  For  a  >  ac  there  is  a  finite  interval  <  x*  <  x*  in  which  the  flow  can 
support  linearly  unstable  Gortler  modes  and  this  interval  increases  with  increasing  a. 

3.  Evolution  equations  for  large  amplitude  Gortler  vortices 

We  now  proceed  to  develop  an  asymptotic  solution  of  the  governing  equations  (2.8a- 
d)  valid  in  the  short  wavelength  limit.  We  assume  that  the  solution  (u(r,  £),  v(x ,  £))  of  the 
classical  boundary  layer  equations  (2.5a,b)  becomes  linearly  unstable  to  Gortler  vortices 
at  x  =  x*  where  x*  is  to  be  determined  from  (2.12). 

In  the  0(h~ */2)  boundary  layer  the  effective  Gortler  number  G  becomes  0(h:i/2)  so 
that,  on  the  basis  of  the  nonlinear  theory  of  Hall  and  Lakin  (1988),  we  anticipate  vortices 
with  wavenumber  k  =  0(h3/s)  (the  scaling  G  =  0(k4)  is  that  which  is  appropriate  to  the 
right-hand-branch  of  the  linear  neutral  curve,  see  Hall  (1982)).  At  first  sight  these  scalings 
appear  to  be  in  disagreement  with  those  presented  in  Hall  and  Bennett  (1986);  in  fact, 
there  is  no  discrepancy  when  one  considers  the  boundary  layer  scalings  (2.7).  The  analysis 
of  Hall  and  Lakin  demonstrates  that  in  the  presence  of  the  nonlinear  vortex  state  the 
flow  regime  breaks  up  into  three  regions;  a  central  core  region,  in  which  the  mean  flow  is 
induced  by  the  vortex  velocity  field,  which  is  bounded  above  and  below  by  thin  shear  layers 
in  which  the  vortex  amplitude  is  reduced  to  zero.  Outside  the  region  of  vortex  activity  the 
mean  flow  is  governed  by  the  unperturbed  boundary  layer  equations  (see  figure  4).  We  now 
proceed  to  develop  the  corresponding  structure  for  the  flow  described  in  §2;  we  emphasize 
that  the  nonlinear  vortex  state  will  be  confined  wholly  within  the  0{h~ 1//2)  boundary 
layer.  Hence,  we  are  anticipating  that  the  only  effect  on  the  onset  of  separation  due  to  the 
presence  of  the  nonlinear  vortex  state  is  that  the  position  at  which  separation  occurs  will 
be  changed;  the  asymptotic  analysis  of  Smith  and  Daniels  concerning  the  removal  of  the 
Goldstein  singularity  still  remains  valid.  As  the  analysis  mirrors  that  of  Hall  and  Lakin  we 
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will  give  only  the  salient  features  of  that  calculation  (the  reader  is  referred  to  that  paper 
for  full  details). 

To  proceed  we  first  consider  the  core  region,  in  which  the  vortex  amplitude  is  chosen 
so  as  to  have  an  0(1)  effect  on  the  mean  flow  in  this  region.  The  appropriate  expansions 
in  the  core  region  are  then  given  by 

u  =  uo  +  h~ 3,/8  {  Uo\  cos  h3^z  +  •••}, 

,  1  ,  J  (3.1) 

v  =  v o  +  h3^8  |  Voi  cos  h3t*z  +  •••>, 

together  with  similar  expansions  for  w  and  p;  here  •  •  •  denotes  terms  which  do  not  enter  into 
the  subsequent  analysis.  Substituting  the  expansions  (3.1)  into  the  governing  equations 
gives,  for  the  vortex  components 

U0i  +  t/ofVoi  =  0?  Voi  +  2uoFxxUoi  =  0,  (3.2a,  6 ) 


which  for  consistency  requires 

2FxxuoUog  —  1. 

Integrating  (3.3)  yields,  for  the  vortex  induced  mean  flow 


(3.3) 


A(x)  +  t, 
uo  =  \l  — - - , 


(3.4a) 


while  from  the  continuity  equation  we  have 


Xu  \  „  }A(x)+t  ,  Fxxx  {A(x)  +  ^3/2 

v0  =  -B(x)-Ax\l — - - -|-  ' 


(3.46) 


here  A(i)  and  B(x)  are  to  be  determined.  Thus,  in  the  core  region  the  mean  flow  is  found 
from  a  solvability  condition  on  the  vortex  amplitude  equations.  At  this  stage  the  vortex 
amplitude  Foi  remains  undetermined;  this  is  found  from  the  mean  flow  component  of  the 
first  momentum  equation.  Thus 


duo  _  duo  d2uo 
uo  — - V  v0 


dx 


di 


-  FF  —  -  — 
1  2di 


duQ 


Vo\ 
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which  making  use  of  (3.4)  can  be  integrated  to  give 


C(x)  + 


V2 

yoi 


=  -Bs 


'A  +  e  Fixx(A  +  zy 


6F2 

XX 


2  y/Fxx(A  +  0 


-FFX£.  (3.5) 


If  we  let  £j(x)  (j  =  1,2)  denote  the  outer  limits  of  the  region  of  vortex  activity  we  then 
have  from  (3.5),  upon  eliminating  C(x),  a  relation  between  £i  and  £2,  namely 


-B 


Fxxx{A+b)2 

GF2X 

Fxzx(A  +  h)2 

sf2ix 


1 

2y/Fxx(A  4  £1 ) 
1 

2  \J  FXX(A  4  £2) 


-  ff,(,  =  -s^Fk 

-  FF,i2. 


(3.6) 


Thus,  the  vortex  amplitude  in  the  core  region  is  reduced  to  zero  as  £  — >  £12!  however  we 
require  thin  shear  layers  centred  on  £12  in  order  to  smooth  out  the  algebraically  decaying 
vortex  component.  We  note  that  at  this  point  of  the  analysis  we  are  unable  to  determine 
the  functions  A(x),  B(x),  £i(.r)  and  £2(:r). 

We  now  consider  the  solution  in  the  shear  layers.  As  noted  by  Hall  and  Lakin,  the 
thickness  of  these  layers  is  found  by  balancing  diffusion  across  the  layer  with  convection 
in  the  stream-wise  direction.  In  the  case  under  consideration  such  a  balance  shows  the 
thickness  of  these  layers  to  be  0(/i-1/4).  We  then  write 


V  =  A'/4«-6), 


and  in  the  shear  layer  we  replace  d/dx  and  ()/d£  by  d/dx  —  hl/4£'2d/dr)  and  hl/4d/dT), 
respectively.  We  will  subsequently  restrict  our  attention  to  the  shear  layer  centred  on  £2; 
an  identical  analysis  holds  in  the  layer  centered  on  £  =  £1. 

From  the  core  solutions  (3.4-5)  the  appropriate  expansions  in  this  shear  layer  are 
found  to  be 

u  =  uq  4  h~l/4Ui  4  h~^2 u 2  4  •  •  •  4 -h~x!2  j(i£oi  4  r1/f4uu  4  •  •  •)  cos  h3^8z  4  •  •  ■  j  , 
v  =  Vo  T  h  1  ^v\  4  •  •  •  4  hl  | (uqi  T  h  1  ^V\  1  4  •  ■  ■)  cos  h3^8  z  4  •  •  •  j  , 
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together  with  similar  expansions  for  w  and  p.  Substitution  of  these  expansions  into  the 
governing  equations  we  find  that,  after  some  trivial  manipulations,  the  leading  order  mean 
flow  terms  are  given  by 

V 


uo  =  —B(x)  -  Ax 

where  we  have  made  use  of  matching  with  the  core  solutions  in  the  limit  tj  >  —  oo.  The 
governing  equations  for  the  leading  order  vortex  amplitudes  (t«oi ,  uoi )  are  found  to  be 
consistent  whereas  the  equations  for  (un,i>ii)  are  consistent  provided  that 


d2v0\  S(x)rjv  01 

drf  3 


+  fi(x)vou 


which  gives  the  governing  equation  for  the  vortex  amplitude  foi  •  Here 


and  /i  is  a  funcion  of  x,  the  precise  nature  of  which  is  not  required  here.  This  equation  is 
a  modified  form  of  the  second  Painleve  transcendent  and  has  been  shown  to  have  a  unique 
solution,  Hastings  and  McLeod  (1980),  with  the  asymptotic  behaviour 


Vo,  ~  Sr/,  tj  — >  -oo,  foi  V  co¬ 


lt  then  follows  that  the  vortex  is  constrained  to  lie  below  f  =  £2  (a  similar  analysis  of  the 
shear  layer  about  shows  that  the  vortex  is  constrained  to  lie  above  £  =  £i)  so  that  the 
nonlinear  vortex  is  confined  to  the  0(1)  core  region.  We  note  from  (3.7a— c)  that  the  mean 
flow  in  the  shear  layer  is  essentially  unaltered  by  the  presence  of  the  vortex  (in  fact,  the 
first  two  terms  in  the  stream-wise  mean  flow  expansion  in  this  layer  are  obtained  from 
the  Taylor  series  expansion  of  the  core  flow  about  £2)-  Hence  the  mean  flow  outside  the 
core  region  must  have  and  v,  in  the  limits  £  — *•  £,  ,£+,  defined  by  the  core  solution 
evaluated  at  £  =  ^1,2  respectively. 
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Outside  the  region  of  vortex  activity  there  is  only  a  mean  velocity  field  and  we  write 


u  =  u  +  0(h  ’/2),  v  =  v  +  0(h  1/2), 


where  (tt,t>),  defined  in  (0,£i)  and  (£2,oo),  satisfy 

_du  _du  d2u 

Ufc+Vd£  =  d^+FFx' 

This  system  is  to  be  solved  subject  to 


du  dv 

fa  +  d(=0' 


u  =  v  =  0,  £  =  0,  u  — y  F,  £  — >  oo, 


(3.9a,  6) 


>1  +  0 


»  ««(£>)  = 


2y/Fxx(A  +  ^y 


v((j)  =  -B(x)-AxJ*±k+F?f. 


0  =  1.2) 


(3.9c  -  e) 


The  system  (3.8,9)  together  with  (3.6)  then  specifies  a  free-boundary  problem  for  £i(x), 
£2(x)  and  the  functions  A(x),  B(x)  together  with  the  skin  friction  u^(x,  0).  The  solution 
of  this  system  is  a  full  numerical  problem  (we  note  that  the  similarity  reduction  available 
to  Hall  and  Lakin  is  not  applicable  here).  Before  proceeding  with  a  discussion  of  the 
marching  procedure  used  to  solve  this  system  we  note  that,  in  order  to  solve  (3.8,9),  it  is 
necessary  to  specify  an  ‘initial  state’  for  the  mean  velocity  field  from  which  the  solution 
can  be  obtained  in  a  systematic  way  by  marching  downstream.  This  requires  finding  an 
asymptotic  form  for  the  solution  close  to  the  position  x*  at  which  the  original  boundary 
layer  becomes  unstable. 

We  assume  that  the  unperturbed  boundary  layer  flow  first  becomes  unstable  at  the 
position  x  =  x*  and  that  the  instability  originates  in  a  thin  layer  centred  on  £  =  £*.  We 
define  I  =  x-x*<l  together  with  the  similarity  variable  rj  =  (£  —  ^*)/X^2  and  expand 
and  £2  in  the  form 


(,=C-  +  o(X), 

(2  =  C  +x''2(  +  o(x), 


(3.10) 
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(these  expansions  follow  from  the  weakly  nonlinear  theory  of  Hall  (1982b)).  Consider  first 
the  core  region;  following  Hall  and  Lakin  (1988)  we  express  (u0,vo)  in  the  form 

u0  =  U  +  X3^2Um(t]) -\ - ,  u0  =  V  +  XVm{t])  H - ,  (3.11a,  b) 


where  (u,  v)  are  the  boundary  layer  velocity  components  in  the  absence  of  the  vortex 
motion.  In  a  small  neighbourhood  of  the  neutral  location  (x*,£*)  we  have  the  Taylor 
series  expansions  for  (u,v)  and  Fxx  given  by 

u  =  u00  +  X'^TjUio  +  X(tj2U20  +  u0i)  4-  X3/2r;un  -f - , 

v  =  v00  +  Xl/2r)v io  +  X(r]2v2o  +  woi)  +  X3/2r)vu  -\ - ,  (3.12a  -  c) 

Fxx  =  Fq  +  XF\  +  •  •  • . 


Substituting  (3.11a,b)  together  with  (3.12a-c)  into  (3.3)  and  equating  like  powers  of  X 


gives 


2-FoUooWio  =  1,  «i0  +  2iioou2o  =  0, 

Pm  =  Xo<i-W,  vm  =  -^-  +  v0m, 


(3.13a  -  c) 


where 


A0  =  —  Uji 


UlQUpi  _  FiUjq 

noo  Fq  ’ 


uiou2o 


and  Vjb  is  a  constant  of  integration  (required  to  ensure  continuity  of  Vm  across  the  shear 
layers).  We  note  that  (3.13a)  is  automatically  satisfied  as  x*  is  the  neutral  position,  while 
(3.13b)  is  satisfied  since  £*  is  chosen  so  that  (2.13)  holds. 

Following  Hall  and  Lakin  (1988)  we  can  show  that  the  vortex  amplitude  V01  vanishes 
at  T)  —  where 

■  12  Co  +  2Ao  /o  1  a\ 

^  “  6Aj  +  u00A0  ’  ^  ^ 


with  Co  an  unknown  constant  to  be  determined.  It  then  follows  from  (3.10)  and  (3.14) 


V-  =  “6  »7+  =  t 

To  determine  r)±  it  remains  to  determine  the  constant  Co. 
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Outside  the  region  of  vortex  activity  we  now  write 


u  =  u  +  X3/2Um(t])  H - ,  v  =  v  +  XVm(v)  4 - ■ 


We  note  that  the  functions  U m  and  Vm  are,  respectively  odd  and  even  functions  of  r\  and 
so  we  will  restrict  our  attention  to  the  region  77  >  tj+.  Here  we  find  that  U m,  Vm  satisfy! 


d?UM  woo 
drj2  2 


dVM  ,  3rr  t]  dU m 

dr]  +2  M  2  dr] 


which  are  to  be  solved  subject  to  the  matching  conditions 


(3.15) 


Um  —  ^0 f?+  —  U'm  —  Ap  —  3Air/^, 


(3.16a) 


and 

Um  — ►  0  t]  — >  00. 

The  appropriate  solutions  are  found  to  be* 

_  (Ao?7+  -  Air^jexp  {-u00(772  -  ?7+)/8}  #(7/2,  y/^rj) 


Um 


U(7/2,^fr,+) 


V 

-2I 


Vm  =  ^ Um  -2  /  UMdr], 


(3.166) 


(3.18a,  6) 


where  U(a,x)  is  the  parabolic  cylinder  function.  From  (3.17a)  and  (3.18a)  we  then  have 
the  eigen-relation^  for  r/+ 


uoo  f  467(9/2,  y/^-T/+)  juoo  )  _  Ap  —  3Ai?7^. 
T  }  P{7/2,  7/WU)  +  V  ~2~'1+]  ~  A„r,+  -  An>  ’ 


(3.19) 


which  then  determines  the  constant  Co  in  (3.14). 

With  the  above  asymptotic  solution,  valid  near  a:*,  we  can  construct  the  composite 
velocity  field  together  with  initial  values  for  the  function  A(x),  B(x),  £1  (x)  and  &(x)  which 
will  then  be  used  to  start  the  marching  procedure  to  be  described  in  the  next  section. 

^Note  the  typographical  error  in  the  corresponding  expressions  of  Hall  and  Lakin  (1988). 

bbid. 
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4.  Numerical  solution  of  the  free  boundary  problem 

The  numerical  scheme  used  to  solve  the  free  boundary  problem  (3.8-9)  is  based  on  that 
used  by  Hall  and  Lakin  (1988);  for  this  reason  we  only  give  the  relevant  details  here.  We 
assume  that  a  solution  is  known  for  x*  <  x  <  x  (which  in  practice  must  be  constructed 
from  the  asymptotic  solution  for  the  initial  development  of  the  nonlinear  vortex  state 
presented  above).  The  scheme  then  advances  this  ‘initial’  solution  downstream  to  x  +  e  for 
some  sufficiently  small  value  of  e. 

In  the  course  of  our  calculations  it  was  found  appropriate  to  solve  a  modified  system 
of  the  form  outlined  in  §2.  Thus  we  define 


u  =  Fu,  v  =  v  —  Fx£. 


The  numerical  scheme  then  proceeds  by  defining  a  new  variable  ip  by 


V>  =  f;  j  =  l,  0<£<6,  j  =  2,  6<£<oo. 

The  modified  form  of  the  system  (3.12)  is  now  written  as 


du  1  d2u 


Fu~  - 


dx  dip 2 


F*(l-u2)- 


( v  -  Fx£j<p)  du  Ffj  du 


+ 


ipu- 


ti  dip  "  dip  ’ 


dv  _  ,  (p  (  ,  Fdu  rpFi'j  du 


dx  dip  J  ’ 


(4.1a,  b) 


where  j  =  1  for  0  <  ip  <  1  and  j  =  2  for  1  <  xp  <  oo.  This  system  is  to  be  solved  on  (0,  oo) 
in  terms  of  ip  with  boundary  conditions 


u  =  v  =  0,  ip  =  0,  u  — >  1  ip  — »  oo, 


(4.2) 


together  with  the  ‘discontinuous’  conditions  at  ip  =  1 


u±  =  -s 


I  lA  +  ti 

F 


ut)±  — 


2  Fy/Fxx(A  +  (jY 


V± 


=  -B  -  A, 


lA  +  Zj  FIXX 

n  I  r> 


(A ±(j 

{  Fxx 


3/2 


+  Fx£j. 


(4.3a, b ) 


15 


Here  the  ±  signs  corresponds  to  the  limits  xp  ►  1+,  0  — ►  1“  while  the  index  j  =  1  is 
associated  with  the  ’  sign  and  j  =  2  with  the  *+’  sign.  The  system  is  then  closed  by  the 
condition  (3.6). 

The  solution  of  (4.1a)  is  advanced  for  xp  in  the  range  0  <  xp  <  1  by  using  the  scheme 

F Xlniin  2  2  (^n+1  2un  -f"  Un- 1)  =  Fun  C.F r(l  XI n ) 

■  (4.4) 

2(  h  ^ Vn  ~~  ~  Fxpn^jXln)  (un+l  —  un  —  l)i 

where  h  is  the  vertical  grid  spacing,  a  tilde  denotes  that  the  quantity  is  evaluated  at 
x  +  e  while  an  index  n  indicates  a  quantity  evaluated  at  the  grid  point  xp  =  nh.  In  order 
to  solve  the  tridiagonal  system  (4.4)  we  must  make  an  initial  guess  for  and  then  set 
=  (£j  —  £j)/e.  In  solving  (4.4)  we  ensure  that  the  boundary  condition  on  ti  at  xp  =  0 
and  (4.3b),  with  j  —  1,  are  satisfied.  The  solution  of  (4.4)  then  yields  un.  We  discretize 
the  modified  continuity  equation  (4.1b)  as 


vn+ 1  Vn  — 

2 h 


Fx(  1  —  un)  + 


F(un  Un)  Fxpn^j  (un-f.]  Un-\) 


which,  given  un  from  (4.4),  allows  us  to  determine  vn  at  the  position  x  +  e  for  0  <  xp  <  1. 
As  the  equation  for  v  is  first  order  in  xp  only  the  boundary  condition  on  v  at  xp  =  0  can 
be  satisfied.  Thus,  although  the  solution  of  (4.1),  for  xp  in  the  range  0  <  xp  <  1,  has  been 
obtained  at  x  +  e  the  additional  boundary  conditions  (4.3a,c),  with  j  —  1,  have  not  been 
satisfied.  These  conditions  can  now  be  used  to  obtain  a  value  for  A  and  an  improved  value 
for  £i  by  writing  them  in  the  form 


li  =  F2Fxzu2_  -  A, 

«  V  Fxx  3  l  Fxx 


+  Fx£  i, 


(4.6a,  b ) 


Here  B  is  the  current  guess  of  B  at  x  +  e.  This  scheme  can  be  applied  in  a  similar  manner 
to  the  region  xp  >  1.  The  equation  (4.4)  is  now  solved  for  u  subject  to  u^>lasxp—>oc 
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together  with  (4.3b)  for  j  =  2,  and  equation  (4.5)  is  solved  for  v  subject  to  (4.3c)  with 
j  =  2.  Writing  (4.3a)  and  (3.6)  in  the  form 

6  =  F2Fxtu\  -  A, 

-  +  ffx{£\  -  &), 

(4.6c,  d ) 

the  first  of  which  then  determines  a  new  value  for  £2  while  the  second  equation  yields  a 
value  for  B.  We  then  have  values  of  u,  v,  A ,  J3,£i  and  £2  at  x  +  e.  This  procedure  is  then 
continued,  making  use  equations  (4.6a-d),  using  Newton  iteration  on  A ,  B,  £  1  and  £ 2  until 
converged  values  of  these  quantities  are  obtained. 

For  definiteness,  our  numerical  calculations  were  restricted  to  the  profile  (2.11).  The 
first  stage  in  the  numerical  solution  of  the  free  boundary  problem  (3.8,9)  involves  deter¬ 
mining  the  neutral  position  from  (2.12,13).  This  was  achieved  using  the  numerical  scheme 
outlined  in  §2;  with  the  initial  conditions  u  =  «  =  0ati=0  the  solution  was  marched 
downstream,  at  each  stream-wise  location  £*  was  determined  from  (2.13),  until  the  point  of 
neutral  stability  x*,  determined  from  (2.12),  was  reached.  At  this  point  an  initial  guess  for 
6(®)j  £2(^)1  A(x)  and  2?(x),  for  x  =  x*  +0.001,  was  obtained  using  the  procedure  outlined 
at  the  end  of  §3.  The  initial  condition  required  for  the  velocity  field  at  x,  namely  u(r, £) 
and  v(x,  £),  were  obtained  using  the  scheme  (2.10)  and  these  values  were  then  interpolated, 
using  cubic  spline  interpolation  in  order  to  ensure  continuous  second  derivatives,  to  yield 
u  and  v  on  a  uniform  grid  in  rp  in  the  intervals  [0, 1)  and  (1,  i/’oo)-  A  similar  interpolation 
procedure  was  used  on  the  vortex  induced  mean  flow  corrections  Um,Vm  thus  allowing  the 
composite  velocity  field  to  be  computed.  The  composite  velocity  field  so  obtained  was  then 
used  as  the  initial  conditions  for  the  marching  scheme  described  above.  The  solution  was 
marched  downstream  until  such  point  as  the  limit  of  the  region  of  vortex  activity  tended 
to  zero  (ie  at  the  point  xf  where  £i(5f)  =  £2(^)1  see  figure  4).  At  this  point  the  velocity 
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field  was  then  used  as  the  initial  condition  for  the  scheme  (2.10)  which  was  then  used  to 
march  downstream  until  such  point  at  which  reversed  flow  is  first  encountered. 

5.  Conclusions 

Before  proceeding  with  a  discussion  of  our  results  we  make  some  comments  concerning 
the  stability  characteristics  of  our  numerical  scheme.  Initially  the  normal  step-size  h  was 
fixed  at  h  =  0.01  with  V’oo  =  16  thus  yielding  100  and  1500  normal  grid  points  in  the 
intervals  0  <  ip  <  1  and  1  <  t/;  <  ifioo,  respectively.  With  h  fixed,  various  calculations 
were  performed  with  different  values  of  the  streamwise  step-length  e;  from  e  =  0.001  to 
(  =  0.00025.  For  the  range  of  e  considered  the  position  xp,  at  which  the  vortex  region 
vanishes,  was  found  to  vary  by,  typically,  0.9%  over  the  value  obtained  for  the  smallest 
value  of  e  considered.  Similarly,  the  normal  position  £f  was  found  to  vary  by  0.7%  over 
the  range  of  e  considered.  We  deem  our  numerical  scheme  to  be  stable  and  for  e  =  0.00025 
suitable  accuracy  has  been  attained.  All  of  the  subsequent  results  were  obtained  for  the 
choice  of  parameters  h  =  0.01,  e  —  0.00025  and  V’oo  =  16. 

Consider  figure  5,  in  which  we  present  a  plot  of  the  positions  (r)  and  &(x)  bounding 
the  region  of  vortex  activity  for  the  constriction  amplitudes  a  =  1.0304  and  a  =  1.0923. 
As  anticipated  from  the  results  of  Hall  and  Lakin  (1988)  we  see  that  the  region  of  vortex 
activity  initially  thickens  until  such  a  point  at  which  the  maximum  normal  extent  of  this 
region  is  achieved.  As  x  is  further  increased,  the  region  of  vortex  activity  constricts  until 
the  point  xp  at  which  £i(rf)  =  ^(^f)-  For  x  >  xp  the  flow  can  no  longer  support  the 
large  amplitude  vortex  state  described  above. 

The  effect  of  increasing  the  amplitude  o  of  the  constriction  is  evident  from  a  com¬ 
parison  of  figures  5a,b.  The  minimum  values  of  in  figures  5a, b  are  2.19  x  10~2  and 
1.08  x  10-2  respectively  whereas  the  maximum  values  of  £2  are  1.384  and  1.752  respec¬ 
tively.  Thus  the  effect  of  increasing  a  is  twofold;  firstly  the  streamwise  extent  of  the  vortex 
region  is  increased  with  increasing  a  (a  consequence  of  the  change  in  the  neutral  positions 
as  represented  in  figure  3)  and  secondly,  the  normal  extent  of  the  vortex  region  increases 
with  increasing  0,  through  both  the  increase  in  the  maximum  value  of  £2  and  a  decrease 
in  the  minimum  value  of 
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Turning  now  to  the  question  as  to  what  effect  the  nonlinear  vortex  state,  present 
in  the  interval  (x*,Xp),  has  upon  the  onset  of  separation  we  present,  in  figures  6a-b, 
plots  of  the  skin  friction  in  the  presence  of  the  vortex  (broken  line)  and,  for  comparison, 
the  corresponding  skin  friction  for  the  unperturbed  boundary  layer  (solid  line),  for  the 
representative  cases  o  =  1.0304  and  o  =  1.0923.  Here  the  skin  friction  is  defined  according 
to 

du(r,0)  F(x )  du(x,0) 


r  =  F(x)- 


(5.1) 


d£  6(z)  dtp 

where  xp  is  the  modified  independent  variable  and  u  is  the  modified  dependent  variable  used 
in  the  numerical  scheme  described  above.  From  a  computational  standpoint  we  replace 
(5.1)  by 


r  = 


F(x)  du(x,  0)  F(x) 


{4 u(x,  h)  —  u(x ,  2 h)}  +  0(h2). 


(5.2) 


Zi(x)  dip  2/i£i(r) 

The  effect  on  the  boundary  layer  due  to  the  vortex  is  readily  seen  through  the  dramatic 
increase  in  the  skin  friction  (over  that  for  the  unperturbed  state).  In  the  presence  of 
the  vortex  the  skin  friction  now  develops  two  maxima,  the  first  being  due  to  the  vortex, 
whereas  the  second  is  seen  to  coincide  with  that  for  the  unperturbed  state.  The  ratio  of 
these  maxima  is  seen  to  increase  with  increasing  a;  for  the  case  of  figure  6a  this  is  found 
to  be  73%  while  for  figure  6b  is  93%.  However,  in  both  figures  6a, b  we  see  that  the  vortex 
induced  skin  friction  subsequently  decreases  to  below  the  value  found  for  the  unperturbed 
state  followed  by  an  increase,  due  to  the  presence  of  the  favourable  pressure  gradient,  until 
such  point  as  the  second  maxima  is  attained  (at  which  point  the  difference  between  the 
vortex  induced  state  and  the  unperturbed  state  is  barely  discernible).  The  skin  friction 
now  decreases  until  the  point  of  separation  is  encountered,  the  location  of  which  is  identical 
for  both  the  unperturbed  state  and  the  case  in  which  the  flow  supports  a  nonlinear  vortex. 

We  are  then  given  to  conclude  that  the  presence  of  a  nonlinear  vortex,  in  a  region 
of  concave  curvature  upstream  of  the  position  of  maximum  channel  constriction,  has  no 
discernible  effect  upon  the  position  at  which  separation  is  encountered.  This  result  is 
somewhat  surprising  in  light  of  the  comments  of  Hall  and  Lakin  (1988)  (pp  443)  that  the 
decay  of  the  vortices  does  not  allow  the  boundary  layer  to  return  to  its  undisturbed  state. 
In  the  present  case  it  would  appear,  by  the  position  at  which  separation  is  reached,  that 
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the  boundary  layer  does,  in  fact,  return  to  its  unperturbed  state.  In  figure  7  we  present 
a  plot  of  the  streamwise  velocity,  for  the  case  of  the  unperturbed  boundary  layer  and  the 
case  when  the  vortices  are  present,  at  the  position  xp. 

In  the  analysis  presented  we  have  ignored  two  other  possible  mechanisms  for  linear  in¬ 
stability  and  subsequent  nonlinear  growth,  namely  Tollmien-Schlichting  waves  and  inviscid 
Rayleigh  modes  (see  Tutty  and  Cowley  (1986)  for  a  discussion  of  the  Rayleigh  instability 
in  interactive  boundary  layer  flows).  The  nonlinear  development  of  these  modes,  together 
with  the  possible  importance  of  the  vortex-wave  interaction  mechanisms  described  by  Hall 
and  Smith  (1988),  could  have  a  dramatic  effect  on  the  results  and  conclusions  of  this  work. 

Finally,  for  the  (somewhat  more  realistic)  case  of  moderately  constricted  channels,  for 
which  h  —  0(1),  the  effective  Gortler  number  is  now  an  0(1)  quantity  and  we  would  expect, 
on  the  basis  of  the  linear  theory  of  Hall  (1983)  for  vortices  with  an  0(1)  wavenumber,  that 
the  flow  will  support  linearly  unstable  Gortler  modes.  The  extension  into  the  nonlinear 
regime  then  involves  the  numerical  solution  of  the  full  nonlinear  governing  equations  as 
was  carried  out  by  Hall  (1988)  for  the  external  boundary  layer  problem.  In  the  present 
case  the  governing  equations  show  distinct  similarities  with  those  derived  by  Denier  and 
Hall  (1992)  to  describe  the  nonlinear  evolution  of  the  most  unstable  Gortler  vortex  mode 
the  solution  of  which  resulted  in  a  separated  flow.  In  the  context  of  an  0(1)  constricted 
channel,  the  resolution  of  this  problem  would  require  the  numerical  solution  of  the  full 
nonlinear  governing  equations.  This  is  currently  the  subject  of  ongoing  research. 
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Figure  1.  A  sketch  of  the  flow  configuration  under  consideration. 
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Figure  2.  Plot  of  skin  friction  r  =  u^(a:,0)  versus  x  for  the  profile  (3.3). 
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Figure  6a.  Plot  of  the  skin  friction  r  in  the  presence  of  the  vortex  (broken  line).  For 
comparison  the  skin  friction  in  the  absence  of  the  vortex  is  presented  (solid  line).  The 
•  represent  the  position  at  which  the  vortex  region  has  decayed  to  zero  (ie  x f).  (a) 
a  =  1.0304,  (b)  a  =  1.0923. 
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Figure  6b.  See  figure  6a  for  legend. 
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Figure  7.  Streamwise  velocity  profiles  for  x  =  Xf.  The  vortex  induced  mean  flow  (broken 
line)  and  the  unperturbed  mean  flow  (solid  line). 
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